Observation of the Presuperfluid Regime in a Two-Dimensional Bose Gas 
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In complementary images of coordinate-space and momentum-space density in a trapped 2D Bose 
gas, we observe the emergence of presuperfluid behavior. As phase-space density p increases toward 
degenerate values, we observe a gradual divergence of the compressibility k from the value predicted 
by a bare-atom model, Kba- k/ Kba grows to 1.7 before p reaches the value for which we observe the 
sudden emergence of a spike at p = in momentum space. Momentum-space images are acquired 
by means of a 2D focusing technique. Our data represent the first observation of non-mean-field 
physics in the presuperfluid but degenerate 2D Bose gas. 
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Because of the enhanced role of fluctuations in low- 
dimensional systems [l|, a two-dimensional (2D) Bose gas 
at nonzero temperature does not have long-range phase 
coherence. In a homogeneous system there can be at best 
only a quasicondensate, no true Bose-Einstein condensa- 
tion (BEC). Under the combined effect of interactions 
and quantum degeneracy, however, there is nonethe- 
less a phase transition known as Berezinskii-Kosterlitz- 
Thouless (BKT) associated with the unbinding of vortex 
pairs [2]. Below the critical temperature Tbkt, the sys- 
tem is superfluid. 

Experiments in 2D atomic gases [3|-[s| are usually con- 
ducted in the presence of an inhomogeneous trapping 
potential. In the complete absence of interactions, the 
confining potential can resurrect a traditional BEC [7j, 
but for realistic experimental parameters, interatomic in- 
teractions tend to suppress BEC by smoothing out the 
spatial profile [ll-lsl, Hill of the mean density to the point 
where the sample can be understood as a collection of 
locally uniform spatial regions, each of which is charac- 
terized by a particular local density and thus a partic- 
ular local value of Tbkt- Although these local regions 
may be too small to test in detail the coherence-related 
predictions of BKT theory^qualitative effects have been 
observed in experiment [3|, [6| . 

Our particular interest is in the region just to the warm 
side of Tbkt- In an earlier experiment on bosons trapped 
in a 2D optical lattice, we observed a proliferation of vor- 
tices as we warmed through the discrete-case equivalent 
of Tbkt [1Q| • But in that experiment a great many meso- 
scopic condensates were present, one at each lattice site, 
on both sides of the BKT transition, because they had 
condensed at a Tbec distinct from and well above Tbkt- 
For the continuous case, in contrast, there is no cor- 
responding second transition temperature above Tbkt- 
But if the cooling gas has by Tbkt already become a 
medium that can support vortices, whether bound or not, 
then heuristically we see that it must have continuously 
evolved from a fully fluctuating nondegenerate gas into 
a sort of presuperfluid with suppressed density fluctua- 



tions pj|. Theory [8l. lllMl5j validates this intuition, and 
experiments [5[ have in turn been consistent with predic- 
tions of that theory. Up until now, however, experiments 
have not been directly sensitive to the properties of the 
presuperfluid, T > Tbkt gas. The goal of the present 
work is to provide a minimal-assumption, first empirical 
look at this exotic regime. We emphasize key features 
of this approach: 1. Our line of sight is along the axis 
of tight confinement: we do not need to do a deconvo- 
lution of our images to get the 2D density distribution. 
Steps are taken to minimize systematic errors in density 
measurements. 2. We analyze our in situ images to ex- 
tract the local compressibility, a quantity directly sensi- 
tive to local microscopic physics. 3. We use a 2D focusing 
technique to record high-resolution 2D momentum-space 
images complementary to the coordinate-space images. 
We make corresponding inferences about nonlocal coher- 
ence. 4. We use a simple but robust "bare-atom" model 
to correct the observed density for the presence of a small 
population in excited states in the tight confinement di- 
rection, and to determine a bare-atom compressibility 
with which to compare our observations. 

Experimentally, we create a stack of well-isolated 
quasi-2D layers by superimposing a one-dimensional, 
blue-detuned optical lattice with lattice spacing 3.8 \im 
onto a magnetically trapped, evaporatively cooled cloud 
of Rb-87 atoms. Within each layer, approximately 
6.2 x 10 5 atoms feel a harmonic potential character- 
ized by frequencies (cj r , cj 2 )=27t(10,1400) Hz. The char- 
acteristic dimensionless 2D interaction strength is g = 
VS7r(a s /aho) = 0.093, where a s is the 3D scattering 
length and is the z harmonic-oscillator length ["Lrjj . 
The atoms are allowed to equilibrate in their 2D geom- 
etry before probing occurs. Right before probing the 
resulting coordinate- or momentum-space distribution, 
we apply a microwave pulse together with a transient 
magnetic field to pump atoms in the central layer into 
another hyperfine state, resonant with the probe light. 
The strength of the microwave selection pulse is adjusted 
to keep the peak optical density of the imaged fraction 
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FIG. 1. (a) A trap- focused, momentum- space image, (b) An 
in-situ, coordinate-space image. The corresponding azimuthal 
averages are shown in Figs. 2(g) and 2(c). 



within the linear dynamic range of our imaging. 

Our momentum-space imaging makes use of a focus- 
ing technique which is an extension to 2D of a proce- 
dure developed for imaging ID momentum distributions 
17 . fl8| . It yields a much cleaner separation of mo- 
mentum and coordinate-space distributions than is ob- 
tained in earlier experiments 0, E[. After selecting a 
single layer, we turn off the optical lattice and let the 
layer expand into a purely magnetic trap with frequen- 
cies (co> r , co^)=27r(5.2,10.4) Hz. Because of the 140:1 as- 
pect ratio of the cloud, the expansion is initially purely 
axial, very rapidly reducing the 3D density while not af- 
fecting the in-plane coordinate- or momentum-space dis- 
tributions. After this near-instantaneous suppression of 
the repulsive atom- atom interactions, each atom under- 
goes free harmonic motion in the x — y plane. After a 
dwell time t = \ just as the initial 2D momentum- and 
coordinate-space distributions have swapped, we take an 
absorption image [Fig. 1(a)]. Scaling the spatial coor- 
dinate by muj r yields the momentum-space distribution 
that existed just as the probe sequence began. We take 
azimuthal averages of the absorption images before plot- 
ting and fitting the data. 



For coordinate-space images the procedure is similar, 
but after the lattice turn-off, we wait only 1 ms before 
taking the absorption image [Fig. 1(b)]. The 2D density 
remains essentially frozen while the 3D density - and 
related imaging artifacts [5[ - are much reduced. 

To extract a signal for many-body physics from our 
data, we compare our data to a fully fluctuating, bare- 
atom model, in essence the Paris group's mean-field, 
Hartree-Fock, local-density model [9]. The mean occu- 
pation of a single-atom state k is given by the Bose- 
Einstein distribution, Nk = e (E k -^/k B T_ 1 ? where Ek is 
the energy of the state. For our system, ksT ^> huj r , 
but ksT ~ hu z . We treat the atomic motion semiclassi- 
cally in the in-plane direction, while preserving discrete 
harmonic-oscillator levels in the z direction. The 2D 



coordinate-space density in the jth axial level is 
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where the free particle dispersion e(p) = J^. The local 
chemical potential for the jth level is given by 



fij(r) = global 
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whereas the intralevel interaction energy is 
0j(r) = 2{Anh 2 a s / m) f jjTij . 



(2) 



(3) 



The relevant mean-field interaction energies depend on 
fji which are the normalized density overlap integrals 
over the axial dimension between densities associated 
with axial quantum states j and I. Interaction energies 
are comfortably less than the axial spacing Hlj z , justify- 
ing our treating the axial wave functions as frozen. We 
define a quantity uqq = (47rh 2 a s /m)foo, such that we can 
write 0o(r) = 2^oo^o = 2(h 2 g/m)no(f). Evaluating the 
integral in Eq. (pQ), we get 



rij (r) = — ln(l — e 
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where the de Broglie wavelength X^b — /niksT '. 
For any given value of ytp lohal and T, nj are determined 
self-consistently. For ksT < Huj z , the model converges in 
just a few iterations. 

The bare-atom model is a no-condensate model from 
which all the many-body effects associated with degener- 
ate bosons has been intentionally omitted: the additional 
factor of 2 in front of the parentheses in Eq. ((2]) and (|3j) 
arises from an implicit assumption that the second-order 
correlation function at zero distance is 2, as it would be 
for fully fluctuating, nondegenerate ideal bosons, and not 
1, as for a 3D BEC. Furthermore, e(jp) = p 2 /2m is the 
dispersion relation for independent atoms moving in a 
mean-field potential. There are no collective excitations 
such as phonons. 

All the same, the bare-atom model should do very well 
where phase-space density pj = rij\ 2 dh < 1, true for j > 
in our system. As for the calculated value of no(r), 
this will begin to fail for po ^ 3.5, but a comparison 
observations with the naive, bare-atom no will allow us 
to quantify the telltale discrepancy. 

An analysis of a coordinate-space image proceeds as 
follows. We measure the integrated density in the z- 
direction with contributions from all populated excited 
axial levels. n rneas (f) — X^ n j(r)> but the interesting 
physics lies in rto(r). We compare the results of the 
bare-atom model to observed n meas and fit the param- 
eters T and pP lohal to the low phase-space (hence well- 
understood) regions of the cloud. With the chemical po- 
tential and temperature obtained from the fit, we can use 
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FIG. 2. (a)-(d) Coordinate-space distributions and (e)-(h) 
corresponding momentum distributions. Two distributions 
in the same row are taken under near-identical conditions. 
The thin black curves give the azimuthal averages of n meas 
from the raw images. The thick blue curves in the coordinate- 
space distributions are the ground-state distribution no after 
correcting for nk>o. The spike in momentum space that first 
appears in (g) has no corresponding dramatic change in coor- 
dinate space (c). The vertical dotted line in (c) represents the 
inverse of the momentum resolution limit indicated in (g) and 
is thus a lower limit on the spatial extent of the coherence of 
the population of low-p (high-coherence) atoms represented 
by the area (about 1.4 % of total) under the spike in (g). 



Eq. (J2J) to evaluate the excited state distributions. Then 
we numerically find a self-consistent solution to get n& 
(k > 0), with the constraint no = n meas — ^Z k> o n k'i see 
Figs. HaW). 

Once the corrected ground-state distribution no is ex- 
tracted from n mea s, we calculate the scaled (by n 2 ) 
isothermal compressibility k at each imaging pixel 



dno/dfio = (dno/dr)(djjbo/dr) 



(5) 



Although pp lohal is not a quantity we can know with 
great accuracy in a model-independent way, d/io/dr ~ 
—d(^muj 2 ,r 2 )/dr = —muj^r is known quite precisely, as 
the contributions to po arising from the mean field of 
axially excited atoms are small and correctable, dno/dr 
is determined from our coordinate-space images. Equa- 
tion ([5]) then gives k at each discrete radius in an image. 
We plot the result vs the local phase-space density po = 
noX% in Fig. [31 We compare k to the value K\y a that the 
bare-atom model would predict at the same density. For 



an observed value of no, we numerically solve the bare- 
atom prediction no = — ln[l — exp(fio — 0o{no)) /ksT)/ \ 2 dh 
for /io, and determine how no changes for small changes 
in /io, and thus extract (Fig. [3j). 

As a test of the local-density approximation that is 
central to our analysis, we determine k using images from 
two very different classes of samples: clouds with T = 
171 nK and central po of about 9 [Fig. [3(a)], and clouds 
with T = 128 nK and central p of about 30 [Fig. [3(b)]. 
The shape of n{p) is the same and, in particular, the 
value of po for which the extracted value of n becomes 
distinguishable from n^a is in both cases about 4. 

A note on our preferred method of fixing of a, the cali- 
bration scale factor that relates n meas (r) to the observed 
optical density profile: we considered (i) an error-prone 
calculation based on optical absorption cross section and 
(ii) a model-dependent (even unto logical circularity) fit- 
ting of a in the image analysis but settled finally on (iii) 
interleaving our data runs with auxiliary measurements 
of very low T clouds in which the atoms are in a near-pure 
Thomas- Fermi inverted parabola with negligible noncon- 
densed wings. In this limit, we assume p 9lobal =uoono(0) 
and thus fix a. This assumption means that our mea- 
sured values of ft, scaled as in Fig. [31 are constrained to 
saturate to an average of 2.0 at very high po. In essence, 
we get accurate measurements of the sample density and 
temperature in the exotic, intermediate regime of po, by 
assuming prior good understanding of behavior in the 
experimentally well-characterized regimes of low degen- 
eracy, mean- field at high T and of high-coherence, pure 
condensates at low T. 

In Fig. [21 we present side-by-side pairs of coordinate- 
space and momentum-space distribution taken in a se- 
quence such that pairs represent images of clouds with 
very similar temperatures and total atom number, such 
that the preimaging values of po(0) are the same between 
pairs to within 10%. As T decreases, there is no obvi- 
ous sudden change in the coordinate-space distribution 
no(r), [Figs. [2(a)-(d)], while in momentum space [Figs. 
[2(e)-(h)] a central spike suddenly emerges at T =171 nK 
[Fig. [2(g)]. The inverse width of this central spike is a 
measure of the spatial extent of the coherent fraction in 
the highest-density region of the cloud. Our momentum- 
space resolution is such that the presence of a resolution- 
limited peak implies that at least some coherence extends 
over a central disk of radius 4.5 /im, or > 10 Xdb- From 
coordinate-space images taken under the same conditions 
for which the coherence spike first appears in momentum 
space, we determine that it happens when the central 
value of po = 8.0(0.7) (this critical value p c is deter- 
mined from looking at many more pairs of images than 
are presented in Fig. [2]). We emphasize that from the 
coordinate-space distribution alone, the identification of 
a transition temperature would require model-dependent 
analysis of the smoothly varying distribution, while with 
access to both distributions at once, we readily see that 
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FIG. 3. Scaled compressibility k vs phase-space density po. 
(a) Measured k extracted from images of samples at T = 171 
nK, as in the image in Fig. 2(c). Black circles are data aver- 
aged over the values calculated from images of three separate 
clouds. The blue curves are Kba calculated from the bare- 
atom model, (b) Same but with k extracted from images of 
samples at T — 128 nK. The paired vertical dotted lines in- 
dicate the location of the jump in coherence discussed in the 
text. 



a modest change in the central phase-space density of 
< 15% causes the distribution at p = to jump by a 
factor of 3. 

What have our observations told us about the 2D Bose 
gas as it cools towards the BKT transition? We can 
say empirically that as po varies from about 7.2 to 8.7, 
we see a dramatic increase of coherence in a range >4.5 
/im, jumping by a factor of 3. The transition may be 
even sharper, but temporal drifts limit resolution. The 
predicted [14| critical value is p c = ln(3S0/g) = 8.3. 

Our most interesting observation is that in warmer 
gases, for po « 4, well before the sudden onset of co- 
herence, we can already resolve that compressibility is 
above what a bare-atom model of degenerate bosons can 
account for. As po reaches its coherence- jump value, 8.0, 
K/ftba has already increased to 1.7 (Fig. [3]). It is natural 
to associate the increase in k) ' Kba with the gradual chang- 
ing of the interaction energy from its value in a fully fluc- 
tuating gas, 2^oo^o, to its fully condensed value, ^oo^o, 
and to draw a corresponding inference about the zero- 
range second-order correlation function. For g = 0.093 



the corresponding theoretical prediction (from Ref. [lij 
and Sec. 3.3 of Ref. |ll|) for n/ Kba at p c would be 1.59, 
in reasonable agreement with our observed 1.7. As a 
caveat, our observation of an anomalous n/uba estab- 
lishes definitively only the breakdown of the bare-atom 
model, which could be partially due to the violation of 
the other key bare-atom assumption, that in-plane exci- 
tations correspond to individual atoms with kinetic en- 
ergy e(p) = p 2 /2m. In any case, our data provide a 
first definitive observation of non-mea-nfield physics in 
the presuperfluid 2D Bose gas. 

We are very pleased to acknowledge useful conversa- 
tions with Z. Hadzibabic, J. Dalibard, W. Phillips, M. 
Holzmann, C. Chin, A. Imambekov, L.-K. Lim and V. 
Gurarie. This work was supported by NSF and ONR. 

Note added. — Recently, a relevant experimental 
preprint appeared (l9| . 
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